clear; clc;
% clc;
%% Setup
InputParam_0        =   struct('XI',1/82.5/4,'Pir_SS',0.02/4,'KAPPA',0.00/4,...
                               'NU_C',1,'NU_L',1,...
                               'TrProb_ext',0.99,'Prob_ext',0.33,...
                               'TrProb_H',0.99,'Prob_H',0.37,...
                               'Fin_DomShare',0.79,'Fin_AdjCost',1,...
                               'Cons_FracT',0.5,'Cons_FracH',0.6,...
                               'B_Total_SS',0.91,'b_lb',-0.29,...
                               'WGap_ext',0.0,'WGap_H',0.0);

[InputParam,PP,SS]  =   SubFun_SteadyState_OutputParam(InputParam_0);

save('Result_4_Benchmark.mat','InputParam','PP','SS');